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ANALYSIS AND DESIGN OF A CANTILEVER-MOUNTED RESILIENT-PAD 


GAS-LUBRICATED THRUST BEARING 
by Izhak Etsion* 

Lewis Research Center 

SUMMARY 

A thrust bearing consisting of pads mounted on resilient, metallic, cantilever 
beams is described and analyzed. Compliance and stiffness of the bearing assembly 
are discussed, and the effects of bearing design parameters on performance are 
shown. After the general analysis, a design example is presented for a flat sector- 
shaped gas bearing. A special case where zero axial movement of the runner can be 
obtained is pointed out. 


INTRODUCTION 

Many modem gas-bearing applications require tolerance to dirt ingestion and 
high-temperature operating capability. The first requirement was met by the develop- 
ment of bearings with flexible or compliant geometry. These are the foil bearings and 
the compliant surface or compliant moimted bearings. In reference 1 the various de- 
signs are briefly summarized and their main disadvantages pointed out. The foil bear- 
ings are less successftil as thrust bearings than as journal bearings. Compliant bear- 
ings that use elastomeric materials are temperature limited. Also in reference 1 a 
new type of bearing is described. This is a resilient-pad thrust bearing that does not 
contain any elastomers in the bearing assembly. It consists of sector- shaped flat pads 
mounted as 5 rmmetrically on resilient foil annular beams. This combines compliance 
with hi^- temperature operating capability. 
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The bearing analysis in reference 1 was simplified by assuming rectangular pads 
and treating the bearing as a slider with constant velocity. The annular support was 
approximated by rectangular beams. 

In the course of experimental work being done at NASA Lewis Research Center it 
was found that the analysis of reference 1 was oversimplified and the bearing failed to 
function under moderate and high loads. Hence a better analysis is needed. As a first 
step the actual flat- pad sector- shaped gas bearing was analyzed (ref. 2) by taking into 
account both pitch and roll of the pad; and results for load capacity, power loss, and 
center- of- pressure location were obtained. 

It is shown in reference 2 that the load capacity of a flat sector- shaped pad is 
maximized whenever the pad tilt results in a uniformly minimum film thickness along 
its trailing e(%e. For that reason and to get better control of the support characteris- 
tics, the annular beams of reference 1 are replaced in this report by strai^t canti- 
lever beams. The bearing concept analyzed thus consists of a sector- shaped flat pad 
mounted on a rectangular cantilever beam, as shown in figures 1 and 2. 


ANALYSIS 
Cantilever Design 

The sector- shaped pad is attached to the cantilever beam along a line parallel to 
the pad trailing edge and at a distance D from that edge (fig. 2). This line serves as 
a pivot about which the pad is tilted. As is shown in reference 2 a pure pitch about 
some radial line is eqmvalent to combined pitch and roll about a given point. It is also 
shown in reference 2 that the radial pivot line that will maximize the load capacity of 
the bearing is its trailing edge. This is also true for any pivot line parallel to the 
trailing edge. The pitch y about the pivot line, the distance D, and the film thick- 
ness hg along the trailing edge completely describe the pad orientation. (All symbols 
are defined in the appendix. ) 

The load W acting on the pad at the center of pressure causes a deflection of the 
cantilever beam, as shown in figure 2. From the strength of materials (e. g. , ref. 3) 
the formulas for the deflection and angle of the end of a cantilever due to a force W 
and a moment M are 
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Substituting the expression for M into equations (1) yields 
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To function properly, the pad must be tilted so that both 6|^ and 
In the case where the center of pressure is behind the pivot line, 
immediately imposes the demand of 


y will be positive, 

X - D < 0. This 
cp 


i>2(D-x^p) (4) 

Before startup, if the pad is loaded, the deflection characteristic should be such 

that the runner touches the pad only at its trailing edge. This provides the wedge- 

shaped film needed for lift-off. In that case the load acts at a distance x =0, and 

cp 

condition (4) becomes 


I > 2D 


(5) 


Rxmner Axial Displacement 

The axial displacement of the runner as a fimction of load and speed is an impor- 
tant parameter in many applications where tight tolerances are necessary. We assume 
tiiat under zero load and speed the pad surface is in contact with the runner and that the 
beam remains vindeflected. The runner displacement 6^, under a load W at any 
speed is 
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where h , the film thickness along the pivot line, is simply 
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Hence 


ar=«b-yi>-h2 


( 8 ) 


For nonresilient geometries, such as fixed pads or rigidly supported tilting pads, 

= 0; and at any nonzero load, h >0. Hence it is clear from equation (6) that 
6^ = -hp and the runner is displaced from its zero- load, zero- speed position. For 
flexible support, however, 0 and dj. can be controlled by the beam design. 
Rearranging equation (8) gives the beam deflection 


6, = 6„ + yD + ho 


(9) 


Substituting equation (9) into (3a) and dividing by (3b) yields 

^r + ^2 3 j 

— + B = l ^ ^ 


y 


2(^cp - 
1 + — 


I 


( 10 ) 


or after algebraic manipulations 
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where a is a dimensionless parameter defined as 
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Selecting values for D/x and Z/x and keeping condition (5), which is l/x^> 

cp cp cp 
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2D/x , give by equation (11) the corresponding value for a. Hence the dimensionless 
runner axial displacement is found from 



1 


( 12 ) 


In a fixed- geometry thrust bearing the loaded runner at zero speed touches the pad 
along the trailing edge. When lift-off occurs the runner moves a distance h2 back 
from its zero- speed position and against the load direction, and thus, S^/hg = -1. In 
a tilting pad bearing, where the pivots are rigidly stpported, the back move:^eiit of the 
loaded nmner due to its velocity must also accommodate the tilt of the pad, and hence 
^i/^2 bearings having compliant support, this backward movement is com- 
pensated by the support deflection, and ^ “i* Thus, the nondimens ional param- 

eter a in equation (12) indicates the amount of compliance of the bearing, ranging 
from zero for rigidly supported fixed pads to any amoimt of compliance as a becomes 
larger than zero. 

With q; as a compliance parameter, it is instructive to examine equation (11) be- 
fore any selection of dimensions I and D is made, hi practical design, the beam 
portion extending beyond the trailing edge is limited by the space available for the 
bearing housing. Therefore, the dimension (I - D) rather than the total length I 
(fig. 2) has to be kept within limits. Rearranging equation (11) in terms of (Z - D)/x 

and D/x^„ we have 
cp 



Keeping (Z - D)/x constant and differentiating a with respect to the variable 
cp 

D/x clearly shows that a has a minimum when D/x = 1. Also, from equa- 
tion (13) it is seen that, for values of D/x close to unity, a is well approximated 

u 

by 
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and is affected less and less by as (Z - D)/x^p becomes larger than 1. 

Figure 3 presents values of the compliance parameter a; as a function of pivot 
location D/x^ at various values of (Z - D)/x^p, as well as the loci of Z = 2D and 
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1=0, It should be recalled that the condition for self-startv^) given by (5) is ^ > 2D, 
Hence practical values of a are confined to the region between these two loci. As 
will be shown in a design example, the clearance parameter yr /h„, the minimvun 
film thickness hg, and the center- of- pressure location are known at any given 

design point. Hence for any desired 6^ a proper value of a can be selected by equa- 
tion (12). It is also clear from equation (12) that a special design is possible where 
= 0. Such a design might be advantageous in some applications and is unique to 
compliant bearings. For zero axial movement of the runner at the design point we have 
from equation (12) 


O' 


/Z!o IIceY' 

V»2 '•o/ 


In general, higher values of a are preferred in order to achieve more compliance. 
The maximum compliance is limited, however, by the space in the bearing housing and 
by the acceptable axial displacement of the shaft. 


Bearing Axial Stiffness 

Another parameter that is of importance in bearing design is the overall axial 
stiffness K. It can be easily shown that in a compliant mounted multipad bearing, low 
stiffoess is associated with better load sharing among the pads. In case of any mis- 
alinement of the runner, the out- of- balance exciting force on an individual pad is re- 
duced and the bearing stability might be improved. 

Referring to the bearing "softness" as the inverse of its stiffness we have 

^ _ d(6p _ d(6p 
K d(W^) n d(W) 

The axial displacement of the runner 6^, depends on the variation in and the axial 
movement of the pad trailing edge. Substituting equation (8) for 6^. we have 

(14, 

K d(W) d(W) d(W) 

Since the last term on the right side of equation (14), namely d(h 2 )/d(W), is always 
negative, it is the other terms, describing the change in the trailing- edge displace- 
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ment, that determine the magnitude of the stiffness. The hi^er the chajige in that dis- 
placement, the softer the bearing. 

The terms ''softness" and "compliance parameter, " which are used in this re- 
port, describe two different features of the bearing. The compliance parameter is 
related to the net axial displacement of a loaded runner due to change in its speed from 
zero to the design value. The softness, or stiffness, is related to an incremental 
change in the displacement due to a change in the load at a constant speed. 

Denoting the trailing- edge displacement by 6^ we have 

^ = (, 5 ) 

d(W) d(W) d(W) 


Differentiating equations (3a) and (3b) with respect to W and substituting into equa- 
tion (15) yields 
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Again after rearranging in terms of 
tion (16) becomes 


^ - D rather than I 


and dividing by 



equa- 
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(17) 


Maintaining {I - D)/x constant and differentiating equation (17) with respect to 

cp 

D/x yields a maximum for d(6.)/d(W) at D/x =0 and a minimum at D/x = 1. 

cp L Cp cp 

In practice, Z - D must be larger than D (see condition (5)) and can be three or 
four times larger, hi that case it is clear from equation (17) that the effect of D on 
the bearing stiffness is negligible compared with that of ^ - D. The same result is 
concluded from equation (13) in regard to the compliance parameter a (see also 
fig. 3). Hence the accurate location D of the pivot line is not crucial, leaving much 
flexibility for manufacturing tolerances. 

So far nothing has been specified regarding the nature of the lubricant or the 
clearance shape. Therefore, the preceding analysis is valid for both compressible 
and incompressible hydrod 5 mamlc bearings of flat, stepped, crowned, or any other 
pad surface form. Moreover, not only sector- shaped pads, but also rectangular 
sliders can be used as the load- carrying members. 

To avoid unnecessary torsional moments, the radial or transverse location of the 
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cantilever beam should be such that the beam centerline is perpendicular to the pad 
trailing edge and passes through the c^ter of pressure. 


DESIGN EXAMPLE 

As an example of the procedtire described, a six- pad gas- lubricated bearing having 
the following dimensions and operating conditions was chosen: 


Outer radius of pad, r^, cm 5 

Inner radius of pad, r., cm 2. 5 

Angular extent of pad, /3, deg 45 

Ambient pressure, p„, N/m 1x10 

Dynamic viscosity of gas, ju, N-sec/m 1. 86x10 

Shaft speed, co, rpm 34 000 

Total load, W. , N 74 

I 2 XI 

Young's modidus of beam material, E, N/m 2. 1x10 


The first step is to determine tiie design parameters for optimtim operation, which 
means the highest value for the minimum film thickness at the design load. Since hg 
appears in the denominator in the definition of the compressibility number A (see ap- 
pendix), the hipest value for h .2 is achieved when the bearing is operated at the low- 
est possible A. This corresponds to a curve for constant A in the design chart, fig- 
ure 4 (from ref. 2), which gives the design load at its extrema. The values of those 
extrema, along with the corresponding optimum clearance parameters, are plotted 
against A in figure 5. 

hi the present design the unit load W^/np^A is 0. 166, which corresponds to the 

maximum unit load at A = 50. From the optimum A and the design speed co, hg is 

obtained. The optimum value of the clearance parameter is also found from figure 5; 

thus, the beam slope y is determined. Also x is found at the design point from 

cp 

figure 6, which is crossplotted from the information on center- of- pres sure location 
given in reference 2. 

In our example, for A = 50 and co = 34 000 rpm, the minimum film thickness hg 

is 1. 4x10"^ centimeter. From figures 5 and 6, at the design point, yr^/hg = 3. 2 and 

X /r = 0. 2; hence, x =1 centimeter. A value of (I - D)/x = 4 was selected be- 
cp o cp cp 

cause of space limitation, and a value of D/x^p = 0 was chosen to achieve low stiff- 
ness. From equation (13), o; is 2.45; from equation (12), 6 ^/h .2 is 0.6; and from 

-3 4 

equation (3b), the beam moment of inertia I is 7.75x10 cm . A beam width b of 
2 centimeters was selected, resulting in a beam thickness t of 3. 6 millimeters. 
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From reference 2 the radial location of the center of pressure r / r is 0. 767, which 

cp o 

determines the radial location of the beam. 


Operating Characteristics at Various Loads and Speeds 

After designing a bearing for optimum operation at a given speed and load, it is 
necessary to determine its operating characteristics over a range of operating speeds 
and loads. 

For a given value of the intercepts of the dimensionless pad load W/APg^A 
with curves of constant A such as those of figure 4 are determined. For each inter- 
cept a value of yr^/hg is obtained. Values of are obtained from curves such 

as those in figure 6. From equation (3b) the beam slope y is calcxilated and hg is 
obtained. A value of co at each intercept point in figure 4 can then be calculated. At 
each such point the values of (I - D)/x^p and H/x^p are determined from the selected 
dimensions I and D, and the nondimensional parameter a is obtained from equa- 
tion (13). Then values of are calcialated by equation (12), and the runner axial 

displacement 6^, is found. Figure 7 shows a plot of minimum film thickness as a 
function of shaft speed for the present beam design at two total loads. Fig^e 8 shows 
the runner axial displacement as a function of shaft speed at various total loads. 

Before startup, if the pad is loaded, the deflection characteristics are such that 
the runner touches the pad at its trailing edge and yr^/hg is infinite. When the runner 
starts to rotate, lift-off occurs and hg increases with speed. At zero speed the run- 
ner, under the load, is displaced from its zero- load position and follows the deflected 
beam. However, as the speed increases the runner is lifted back towards its zero- 
load position. And as the speed continues to increase, the runner passes the zero- load 
position. 

The bearing stiffness K at any given speed and load can be easily found from fig- 
ure 8. This is done by crossplotting against 6^, for the given speed and finding 
the slope d(W^)/d(6p at the given load. In figure 9 such variation of 6^, with is 
shown for various shaft speeds. The bearing stiffness is constant over a wide range of 
operating conditions. 

In many practical applications the load on the bearing varies with speed. If this 
variation is known, the corresponding points can be plotted in figure 8 to describe the 
runner axial movement. A particular case is the line 6^ = 0. As was mentioned, a 
special design where 6^ is zero at the design point can be performed by choosing a 
proper compliance parameter a. However, if for such a bearing the off- design load 
can be controlled and related to the speed in order to result in that particular locus, 
a unique design can be accomplished where the runner never moves axially from its 
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zero-speed, zero- load position. 

The friction torque F / to at each point is calculated from the information contained 
in figure 10 (from ref. 2). At each intercept point in figure 4 the values of W, A, and 
^^o^^2 known. Hence from figure 10 the corresponding values of F/Wwhg are 
obtained. The values of to and hg that were already calculated are now used to ob- 
tain the power loss F or the friction torque F/co. In figure 11 the torque is plotted 
against shaft speed at two different loads. 


Summary of Beam Design Procedure 

The beam design procedure can be summarized as follows: 

(1) For a given load and a given bearing geometry, determine the pad unit 
loading W/p^A, where W = Wj./n and A = j8 ^r^ - r?^ j 2. 

(2) From curves such as figure 5, determine the optimum A and the optimum 
clearance parameter yr /h«. 

/ 2 / 

(3) Calculate h 2 at the design speed to from hg = 

value to calcvilate the pad tilt y from the optimum clearance parameter. 

(4) From curves such as figure 6 determine x at the design point. 

cp 

(5) Considering space limitation, select i - D and choose D. (A value of 
D/x^p = 0 will result in the lowest stiffness. ) Use equation (11) to determine the com- 
pliance parameter a and calculate the runner axial displacement from 6 = o;yx - 
h 2 . Check if this displacement is acceptable and alter (I - D)/x^p accordingly to ob- 
tain other values for a and 6„. The term (I - D)/x is mostly limited by the avail- 

JL ^P 

able space in the bearing housing and not by the runner displacement. 

(6) With the final values of I and D, use equation (3b) to calculate the beam mo- 
ment of inertia I. Select a value for the beam width b and calculate its thickness t 
from (12 l/b)^/^. 

A flow diagram summarizing the beam design procedure is presented in figure 12. 
Another diagram for off- design calctilation is given in figure 13. 


SUMMARY OF RESULTS 

The procedure for optimizing the design of a cantilever- mounted, resilient- pad, 
gas- lubricated thrust bearing for operation at a given load and speed was shown to be 
straightforward. A compliance parameter was defined that determines die axial dis- 
placement of the runner and indicates the bearing compliance compared with rigid 
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configurations. The operating characteristics of an optimized bearing design were de- 
termined for a range of operating speeds and loads in a design example. And a way to 
calculate bearing stiffness was described. The effect of cantilever beam dimensions 
on bearing characteristics was discussed, and it was concluded that bearing compliancy 
and stiffness are mostly affected by the length of that portion of the beam that extends 
beyond the pad trailing edge. Increasing this length increases the compliance and re- 
duces the stiffness. It was also shown that a special thrust bearing can be designed 
where the runner axial movement is zero, provided off- design loads can be controlled. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 31, 1976, 

505-04. 
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APPENDIX - SYMBOLS 


pad area, i3^r^ - j/s 
beam width 

distance of pivot line from trailing edge (positive if pivot line is ahead of trail- 
ing edge) 

Yoimg s modulus of beam material 
power loss 

film thickness along pivot line 
minimum film thickness 

O 

moment of inertia of flexible beam, bt /12 

bearing axial stiffness 

length of beam 

moment acting on beam end 

number of pads 

ambient pressure 

radial center- of-pressure location 

pad inner radius 

pad outer radius 

thickness of flexible beam 

pad loading, W^n 

total load 


distance of center of pressure from trailing edge 


compliance parameter, ^ + 1 W_o _2E 




^2 ^o 


angular extent of pad 

slope at beam free end, equal to pad tilt about its pivot line 
deflection of flexible beam 


axial displacement of runner, in direction of load 
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axial displacement of pad trailing edge, in direction of load, - yD 
angular center- of- pressure location, measured from pad leading edge 

compressibility number, Sjuair^/p^hg 

dynamic viscosity 
shaft speed 
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gure 3. - Compliance parameter as function of pivot location at various 
(l-D)*/x^p ratios. 


Unit load, W/ApgA 



2^ 4 6 8 10 12 W 16 18 


Clearance parameter, r^Ql^2 

Figure 4. - Load as function of clearance parameter at various compressibility 
numbers. Ratio of inner to outer pad radius, 0.5; angular extent of 
pad, p, 45^, (From ref. 2. ) 



Center-of-pressure location 



Figure 5. - Maximum available unit load and optimum clear- 
ance parameter as functions of compressibility number. 
Ratio of inner to outer pad radius, r^/rQ, 0.5; angular 
extent of pad, p, 45°. 



Clearance parameter, yrQ/h 2 


Figure 6. - Center-of -pressure distance from trailing edge as function of clear- 
ance parameter at various compressibility numbers. Ratio of inner to outer 
pad radius, Tj/r^, 0.5; angular extent of pad, p, 45° 
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Power loss coefficient, F/Wo;h 
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Figure 9. - Runner axial displacement as function of total load for vari- 
ous shaft speeds. Length of beam, I, 4 centimeters; pivot location, 

D, 0; moment of inertia, I, 7.75xl0"^ cm^. 



Clearance parameter, yrQ/h 2 

Figure 10. - Power loss coefficient as function of clearance 
parameter at various compressibility numbers. Ratio of 
Inner to outer pad radius, 0.5; angular extent of 
pad, p, 45° (From ref. 2 . ) 
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